Integrating the effects of driving forces on ecosystem services into ecological management: A case study from Sichuan Province, China

Driving forces are the factors that lead to the observed changes in the quantity and quality of ecosystem services (ESs). The relationship between driving forces and ESs involves considerable scale-related information. Place-based ecological management requires this information to support local sustainable development. Despite the importance of scale in ES research, most studies have only examined the association between ESs and their drivers at a single level, and few studies have examined this relationship at various scales or analyzed spatial heterogeneity. The purpose of this paper is to explore the significance of the scale-dependent effects of drivers on ESs for localized ecological management. The biophysical values of ESs were calculated using several ecological simulation models. The effects of driving forces on ESs were explored using the geographically weighted regression (GWR) model. Variations in the effects of driving forces on ESs were examined at three scales: provincial, ecoregional, and subecoregional scales. Finally, canonical correlation analysis was used to identify the major environmental factors associated with these variations in each ecoregion. Our results show that (1) the distribution of soil conservation and water yield is highly heterogeneous; (2) four driving forces have significant positive and negative impacts on soil conservation and water yield, and their effects on the two services vary spatially (p < 0.05); (3) the impacts of drivers on ESs vary across different spatial scales, with a corresponding shift in the related environmental factors; and (4) in the study area, at the provincial scale, physical, topographical, and biophysical factors were key factors associated with the variations in the relationship between ESs and drivers, and at the ecoregional and subecoregional scales, physical, socioeconomic, topographical, and biophysical factors all contributed to these changes. Our results suggest that significant differences in topographical conditions (e.g., altitude, slope) can be incorporated for exploring the relationship between drivers and ESs and optimizing ecological management at the provincial scale, whereas significant differences in physical and socioeconomic conditions (e.g., urbanization levels, human activity, vegetation coverage) are more meaningful for localized ecological management at the ecoregional and subecological scales. These findings provide a basis for understanding the relationship between drivers and ESs at multiple scales as well as guidelines for improving localized ecological management and achieving sustainable development.

coordinate this relationship to manage the service to produce desired outcomes at the state scale. In this context, understanding the effects of driving forces on ESs from a multiscale perspective, identifying factors related to them, and determining how they are related would help us to optimize scale-appropriate decisions on localized management and strengthen sustainable development [32,33].
There is increasing awareness of the importance of scale effects for ES analysis and valuation of ESs [34]. Scale effects usually refer to the differences in ecological patterns and processes due to the scale of observation and contain two dimensions, spatial and temporal [3,28,35]. Previous studies related to the scale effects of ES research have mainly focused on the following aspects: supply, demand, and measurement [3,10,31,36]. The scale effects of ES supply can be embodied in the spatial pattern of individual ES [3,5,31], ES interactions [37][38][39][40][41], and ecosystem service bundles (ESB) [3,31,42]. The ES demand by stakeholders at different spatial scales, such as municipal, provincial or national, and international levels, has its own emphasis. Specifically, local stakeholders may emphasize local heritage, culture, or facilities, while national and/or global stakeholders may pay greater attention to the conservation of nature and biodiversity [28]. The ES supply-demand relation is even more complex [43], and such a relation at different scales may result in a scale mismatch [19]. In addition, the scales of observation determine the ES output [42]. Specific scales, such as regular grids or administrative districts (e.g., the township scale and county scale), are widely applied as the research unit [3,5,6]. Despite recognizing the importance of scale effects in ES research, most studies have only examined the associations between ESs and their drivers at a single scale, such as regular grids or administrative districts [6], and few have examined these relationships at different scales [10] or analyzed spatial heterogeneity [7,31].
The analysis of scale dependency is especially important to the ecological management of Sichuan Province.
As an important ecological barrier in the upper reaches of the Yangtze River, Sichuan Province faces huge challenges, especially the issue of ecological sustainability [33,44]. Sichuan Province is characterized by significant regional heterogeneity, which is related to the complexity of geomorphic conditions such as mountains, plateaus, hills, basins, and plains and the diversity of climate zones, including the Plateau Climate Zone, North Subtropical Zone, and Middle Subtropical Zone [44]. Previous research has shown that the impacts of drivers on ESs vary widely across different spatial scales in this area. For instance, [45] found that climate and topography were the main factors affecting the spatial distribution of soil conservation at the provincial scale, while vegetation factors played a significant role at the local scale. Hence, to maintain the vital functions of ESs and promote sustainable development of Sichuan Province, it is necessary to examine the relationships between ESs and driving forces at multiple spatial scales and provide references for regional ecological management.
Through a case study of Sichuan Province, China, this study investigates the significance of the scale dependency of drivers on ESs for localized ecological management. The biophysical values of ESs were calculated and mapped in 2018. The geographically weighted regression (GWR) model was utilized to explore the effects of driving forces on ESs. Variations in their effects on ESs were examined at provincial, ecoregional, and subecoregional scales. Finally, canonical correlation analysis was used to identify the major factors associated with these variations at each spatial scale. The specific objectives of this study were to (1) analyze variations in the effects of driving forces on ESs at three spatial scales; (2) identify the main environmental factors associated with these variations at each spatial scale; and (3) discuss the implications of regional ecological management by analyzing the relationships between environmental factors and drivers.

Study area
Sichuan Province is located at 97˚21' E to 108˚31' E and 26˚03' N to 34˚19 N, Southwest China, covering a total area of 4.86 × 10 5 km 2 [33] (Fig 1). The elevation of the study area ranges from 212 m to 6,904 m and gradually decreases from the western to the eastern side of the study area [46]. Sichuan Province is made up of plateaus and mountains with elevations exceeding 3,000 m in the west and basins and hills with elevations ranging from 500 to 2,000 m in the east [44]. There are significant differences in climate characteristics among different regions. The eastern area is characterized by a subtropical humid monsoon climate with an annual temperature of 16-18˚C and annual precipitation of 900-1,200 mm. The western region is characterized by an alpine climate with an annual temperature of 4-12˚C and annual rainfall of 500-900 mm [45]. In 2018, forestland accounted for 45.57% of the total study area, followed by grassland (25.13%), cropland (13.86%), other types of land (7.88%), built-up land (3.93%), water bodies (2.12%) and garden land (1.5%). The spatial distribution of the main land-use types showed significant regional differences [33]. Known as one of the most important industrial centers in western China, Sichuan Province is currently undergoing rapid urbanization [44]. In 2018, the total gross domestic product of this area reached 40,678.13 billion yuan, and the total population was approximately 83.41 million persons, accounting for 4.4% of the total gross domestic product of China and 5.98% of its national population [46].

Data sources
The data used in this study were obtained using the following sources.
(1) Land use/cover (LULC) data in 2018 with a spatial resolution of 30 m were retrieved from the Resource.
and Environment Science and Data Center (https://www.resdc.cn/data.aspx?DATAID= 264). In this study, the land-use types were classified into six types, cropland, forestland, grassland, waterbody, built-up land, and bare land, based on the "Chinese Classification Criteria of land use" (GB/T21010-2007).
(2) Digital elevation model (DEM) data with a resolution of 90 m were obtained from the Resource and Environment Science and Data Center (https://www.resdc.cn/data.aspx? DATAID=284).
(3) The normalized difference vegetation index (NDVI) with a resolution of 1000 m was retrieved from the Resource and Environment Science and Data Center (https://www. resdc.cn/data.aspx?DATAID=257).
(4) Meteorological data from 42 meteorological stations within the study area were collected from the Resource and Environment Science and Data Center (https://www.resdc.cn/data. aspx?DATAID=230).
In this work, we referred to the environmental factors as spatial variables with high spatial heterogeneity. According to the environmental features and data availability of Sichuan Province, ten environmental factors in the economic, social, ecological, and environmental aspects were selected [46]. Table 1 gives a brief description of these environmental factors. The physical and topographical factors were calculated in ArcGIS 10.6 software (Esri, US) [48]. Socioeconomic factors were retrieved from the Resource and Environment Science and Data Center (https://www.resdc.cn/DOI/doi.aspx?DOIid=33 and https://www.resdc.cn/DOI/DOI.aspx? DOIid=32). In addition, the driving forces of ESs, including gross domestic product (GDP), population density (POP), runoff, and NDVI, were also selected from these environmental factors based on previous studies [26,[49][50][51]. The spatial distributions of the four drivers are displayed in Fig 2. Finally, to unify the resolution of differences and improve the accuracy of the results, all the data were interpolated and resampled to 1000 m resolution by using ArcGIS 10.6 software (Esri, US).

Ecosystem service quantification
Choosing the appropriate indicators to quantify ESs is difficult because the main ESs differ from ecosystem to ecosystem. There is a need to select ESs that are most representative of the ecological problems facing the study area [51]. Sichuan Province has experienced long periods of soil erosion [45] and water shortage [52], which result in soil fertility decline, habitat loss, and local land degradation. Therefore, considering the main ecological challenges, strong linkage to the study area, and data availability, two ESs, namely, soil conservation (SC) and water yield (WY), were selected and estimated in this work.
Soil conservation. Soil conservation (SC) is a service that describes the ability of ecosystems to reduce soil erosion and mitigate the sedimentation of rivers, lakes, and reservoirs through the interception, absorption, infiltration, and fixation of root systems [53]. SC can be calculated as the difference between potential soil erosion and actual soil erosion [45,51]. In this study, the Revised Universal Soil Loss Equation (RUSLE) was used to estimate SC [54]. This equation is expressed as: where SC represents the amount of annual SC (t ha -1 ); R is the rainfall erosivity index (MJ mm ha -2 ha -1 ), which is calculated by using Wischmeier's monthly scale formula [55]; K is the soil erodibility index (t h MJ -1 mm -1 ), which is calculated by using an equation provided by the United States Department of Agriculture [56]; LS is the topographic factor (dimensionless); C is the vegetation cover factor (dimensionless, ranging from 0 to 1); and P is the conservation practice factor (dimensionless, ranging from 0 to 1).
where P is the average annual precipitation (mm), P i is the average monthly precipitation (mm), K is the soil erodibility index, and SAN, SIL, CLA, and C are the contents of sand, silt, clay, and organic carbon, respectively. Water yield. Water yield (WY) is one of the most valuable services to society and an integral ecosystem component that regulates living biomass, the carbon cycle, and the energy budget [8]. WY is a measurement of the total amount of water that flows out of a drainage basin within a specified period of time [8] and is widely determined by the difference between precipitation and evapotranspiration [21,39]. In this work, the InVEST (Integrated Valuation of Ecosystem Services and Trade-offs) model, which is widely used for ES evaluation and mapping, was used to assess WY [57]. Its calculation formula is as follows: where Y X is the annual water yield (mm), P x is the annual precipitation (mm m -1 ), AET is the actual evapotranspiration (mm m -1 ), ω x is a dimensionless ratio of plant-accessible water storage to precipitation, R x is the aridity index for certain land use/cover types, K x refers to the vegetation evapotranspiration coefficient associated with land-use/cover types, ET 0 represents the potential evapotranspiration calculated by the algorithm of the FAO [58], z is a hydrogeological constant that is computed by the method of [59] and AWC is the volumetric plant available water content.

Driver effect analysis
In this study, the geographically weighted regression (GWR) model was adopted to identify the spatially nonstationary relationships between ESs and driving forces. The GWR model is a local linear regression method that extends ordinary least squares (OLS) regression by allowing the regression coefficients to vary continuously over the study area and generating a set of local-specific coefficients for each sample site [60][61][62][63]. Unlike the OLS model, which assumes stationarity across the study area and reflects the "average" or "global" condition of the relationship, the output of the GWR model has proven to be an effective visualization tool to understand spatially unstable variables [8,46,51,64]. Given its strong advantage, the GWR model has important implications for decision-making in ES management [19] and has been widely used to account for spatially nonstationary relationships between ESs and driving forces [8,51,[64][65][66][67]. The GWR model can be expressed as follows [63]: where y i is the dependent variable of sample point i; (u i , v i ) refers to the geographical coordinates (latitude and longitude) of sample point i; β 0 (u i , v i ) is the intercept value for sample point i; β k (u i , v i ) denotes the local regression coefficient of independent variable x ik ; and ε i is the random error or residual for sample point i. In the GWR model, the influence of each observation is determined by its geographic location and its distance from the observed points. Following earlier work, the Gaussian distance decay function was selected to express the spatial interactions as the decay function, which could be calculated by the following equation [68].
where w ij is the geographical weight of sample point j; d ij is the distance between sample points i and j; and h is the kernel bandwidth, which controls the effect of the distance on the weight value. The golden selection search function was used to select the optimal bandwidth [69]. In MGWR 2.2.1 software, four driving forces (GDP, POP, runoff, and NDVI) were chosen as independent variables, and the biophysical values of two ESs (soil conservation and water yield) were chosen as dependent variables.

Scale effect analysis
Based on the maneuverability of management policy, the differentiation laws in physical geography, and previous research, three spatial scales were chosen for this study [70]. The spatial distribution of the three spatial scales is shown in Fig 3. Based on relevant research [3,67], the variations in the effects of drivers on ESs were calculated for each scale by the following steps. First, vector map layers were created at the ecoregion (E) and subecoregion (SE) scales using clip tools in ArcGIS 10.6 software (Esri, US). Thus, we were able to extract the statistical values of the regression coefficients of four driving forces (GDP, POP, runoff, and NDVI) from each of the four ecoregions and 12 subecoregions. Following this, the regression coefficients of each driving force within each vector map were then standardized to a fixed scale (for instance, 0-1 in this study) by the min-max normalization technique for convenient statistical analysis. Finally, the mean value and standard deviation of the regression coefficients of the four driving forces in each ecoregion and subecoregion were calculated using the Zonal Statistics tool in ArcGIS 10.6 software (Esri, US).

Analysis of environmental factors
Finally, canonical correlation analysis (CCA) was used to determine the potential correlations between environmental factors and the regression coefficients of driving forces at three spatial scales. As a multivariate statistical analysis, CCA is a standard tool used to determine the consistency (measured by correlation) between two sets of variables [71,72]. There is a rationale behind this method that is based on finding a pair of linear combinations of U 1-1 and V 1-1 (a linear combination of the dependent variables (X 1 , X 2 , . . ., X k1 ) and a linear combination of the independent variables (Y 1 , Y 2 , . . ., Y k2 )) that have the largest correlation between U 1 and V 1 . The process of finding a pair of linear combinations of U 1-2 . . .U 1-d and V 1-2 . . .V 1-d continues until no correlation between U 1 and V 1 remains statistically significant [10]. CCA represents a high-dimensional correlation between the variables (X 1 , X 2 , . . ., X k1 ) and (Y 1 , Y 2 , . . ., Y k2 ) with a few pairs of canonical variables. In this study, the absolute values of canonical loadings from small to large indicate the internal relationships between environmental factors, and the regression coefficients of driving forces range from weak to strong. CCA was performed using SPSS 26.0.

Spatial distribution of ecosystem services
The distribution of soil conservation (SC) and water yield (WY) is uneven in Sichuan Province, and there are significant differences in the distribution between regions. Low values of SC mean that ecosystems are less capable of reducing soil erosion and mitigating sedimentation in rivers, lakes and reservoirs, while high values of SC mean that this capability is greater. The low values of SC are mostly located in the southeastern and northwestern parts of Sichuan, particularly in the Northwest Sichuan Plateau and Central Hill, while the high values are mainly concentrated in northeastern Sichuan, including the Northern Sichuan Basin and Northern Hill.
Low WY values imply a smaller amount of water flowing out of a drainage basin within a specified period of time, while high WY values imply a larger amount of water flowing out of a drainage basin within a specified period of time. WY values are low in northwestern Sichuan (e.g., Northwest Sichuan Plateau, Alpine Plateau), while high values are mainly found in northeastern Sichuan, particularly in the Mountains around the Sichuan Basin and Low mountains and hills. The distribution of WY in Sichuan Province displays a decreasing trend from southeast to northwest, as well as strong geomorphological distribution characteristics.

Spatial variations in the effects of driving forces on ecosystem services
The spatial variations in the effects of driving forces on ESs were identified by mapping the local regression coefficients of each driving force from the GWR model using ArcGIS 10.6 software (Esri, US). Overall, Fig 4 demonstrates that the four driving forces (GDP, POP, runoff, and NDVI) have significant positive and negative impacts on soil conservation and water yield, and their effects on the two services vary spatially. Fig 4A and 4B shows the effects of GDP and POP on ecosystem services. Both stronger positive and negative impacts of GDP and POP on soil conservation and water yield are mainly observed in the Western Sichuan Plateau. Fig 4(C) shows the effects of runoff on ecosystem services. A stronger positive impact of runoff on soil conservation is observed in the southeast of the province, mainly in Central Hill and Northern Hill, while negative effects of runoff on soil conservation and water yield are mainly found in the northeast fringe. Fig 4(D) shows the effects of NDVI on ecosystem services. NDVI has a significant negative effect on soil conservation in the central and western parts of the province and a positive effect on soil conservation and water yield on the east side of the province.

Scale dependency of driving forces
The variations in the statistical values of the regression coefficient of four drivers within each ecoregion or subregion are shown in Tables 2 and 3. Statistical results indicate that the impacts of driving forces (GDP, POP, runoff, and NDVI) on two services (soil conservation and water yield) differ significantly across three spatial scales (provincial, ecoregional, and subecoregional). In general, regional differences may be small within ecoregions or large between them. However, within some ecoregions (e.g., the Western Sichuan Plateau), there are also differences in the driving force coefficients. In the Southwestern Mountain, the driving coefficients of soil conservation vary significantly between subecoregions, while those of water yield are not significantly different.

PLOS ONE
Integrating the effects of driving forces on ecosystem services into ecological management

Relationships between environmental factors and drivers
Based on the absolute values of canonical loadings from CCA, the associations between environmental factors and driver regression coefficients were identified across three spatial scales (Tables 4 and 5). According to previous research [10], we chose a significance level of 0.05 and excluded correlations that were not significant at this level. At the provincial (P) scale, the key factors underlying the regression coefficients of drivers are the relief degree of the land surface, digital elevation model, proportion of grassland, normalized difference vegetation index, and proportion of built-up land.
At the ecoregional (E) scale, the major factors correlating the regression coefficients of drivers are population density, the proportion of forestland, digital elevation model, the relief degree of the land surface, and normalized difference vegetation index in the Western Sichuan Plateau (WSP); the relief degree of the land surface, the proportion of forestland, normalized difference vegetation index, and digital elevation model in the Southwestern Mountain (SM); the proportion of cropland, digital elevation model, the relief degree of the land surface, and gross domestic product in the Mountains around the Sichuan Basin (MSB); and the proportion of cropland, digital elevation model, the relief degree of the land surface, gross domestic product, and the proportion of built-up land in the Low mountains and Hills (LMH).
At the subecoregional (SE) scale, the primary factors underlying the regression coefficients of drivers are the normalized difference vegetation index, the proportion of forestland, proportion of built-up land, the relief degree of the land surface, and gross domestic product in the Table 4

Physical factor
Socioeconomic factors    Briefly, as shown in Tables 4-6, it is evident that the major factors associated with the regression coefficients of drivers are significantly different at the three spatial scales (provincial, ecoregional, and subecoregional) in Sichuan Province. Specifically, physical, topographical, and biophysical factors were key factors associated with the variations in the relationship between ESs and drivers at the provincial scale, and physical, socioeconomic, topographical, and biophysical factors all contributed to these changes at the ecoregional and subecoregional scales.

Scale dependency of driving forces
There has been some evidence that the impacts of drivers on ESs may vary with scale [6,7,20]. Based on Tables 2 and 3, the scale dependency of the impacts of driving forces on the two services was confirmed in this study. The reason for this shift is that landscape multifunctionality is more evident in larger spatial units than in smaller ones, facilitating a shift in the relationship between services from trade-offs to synergies [3,31,42]. Xu et al. 2017 [5] found that the ES of food supply and tourism resulted in trade-offs at fine spatial scales (<7 km); however, such trade-offs became synergies at a coarser spatial scale (>7 km). Our research also found the same phenomenon between soil conservation and water yield, where this relationship changed from trade-off (at the SE scale) to synergy (at the P and E scales) ( Table 6). The study supports the hypothesis that an increase in the spatial scale of observation will result in a homogenization of the landscape [42].

Relationships between environmental factors and drivers
Sichuan Province. The results in Table 4 reveal that at the provincial (P) scale, the digital elevation model and the relief degree of the land surface were the key factors associated with the variations in the regression coefficients of the four driving forces in Sichuan Province. Topography, as a fundamental physical element, plays an essential role in determining the distribution of land use [8,73,74], human activities [26,75], vegetation [76][77][78], and water resources [79] across the world. In Sichuan Province, there are a plethora of landform types, including plateaus, mountains, hills, basins, and plains. In the east, there are basins and hills, with elevations mostly between 1000 and 3000 m. In the west, there are high mountain plateaus and mountainous landscapes, with elevations mostly above 4000 m. Huge differences in topographical conditions first cause variations in environmental aspects (e.g., temperature, precipitation, wind speed) and then result in remarkable regional social-ecological differences Table 6. The relationship between soil conservation (SC) and water yield (WY) at three spatial levels. "+" represents synergy, and "−" signifies a trade-off.

Scales
Relationship  [44]. As shown in Fig 2, the four drivers in Sichuan Province are extremely unevenly distributed. For example, the highest mean values of population density and gross domestic product are mostly located in the Sichuan Basin, and the lowest values are mainly concentrated in the Northwest Sichuan Plateau ( Table 7). The spatial distribution of the main LULC types exhibited regional differences in Sichuan Province [33]. Grassland and woodland are mainly distributed in western areas, while cropland and construction land areas are mainly concentrated in the eastern and southeastern regions. In addition, the spatial distribution of vegetation is also closely associated with topography: farmland and human settlements are usually found on flat terrain, while forested areas occupy steeper slopes in Sichuan Province. Hence, these results suggest that significant differences in topographic conditions can be taken into account when exploring the influences of drivers on ESs and optimizing ecological management in Sichuan Province. Western Sichuan Plateau. Based on the results of Tables 5 and 6, the proportion of forestland, the proportion of grassland, and the normalized difference vegetation index are the main factors responsible for the variation in the regression coefficients of the four driving forces in the Western Sichuan Plateau. The reason may lie in the fact that grassland and woodland mostly cover western Sichuan Province (especially the Tibetan Autonomous Prefecture of Garz, Tibetan and Qiang Autonomous Prefecture of Ngawa, and Yi Autonomous Prefecture of Liangshan) and have important economic and ecological values [33]. On the one hand, woodland and grassland can provide humans with timber, energy, food, and other products and are the basis for local economic development and population agglomeration. On the other hand, natural ecosystems play a key role in local ecological conservation. For instance, surface runoff can be reduced due to the retention effect of the forest canopy and litter [38]. Therefore, significant differences in the economic and ecological values of local natural vegetation can be considered in managing the relationship between driving forces and ESs and in improving the localized landscape and ecosystem management in the Western Sichuan Plateau. Southwestern Mountain. Based on Table 4, it is evident that the proportion of forestland, digital elevation model, and the relief degree of the land surface account for the majority of the variation in the regression coefficients of the four forces in the Southwestern Mountain. Southwest Sichuan is mainly covered by woodland and characterized by continuous karst topography, sheet erosion, and gully erosion [50,80]. Furthermore, West Sichuan is a region prone to natural hazards [44]. For example, the Central Mountain and Liangshan Plateau are areas with complex geological conditions, highly erodible water, and frequent natural disasters. Heavy precipitation and thin soil layers will increase the frequency of soil erosion, landslides, debris flows, and other natural hazards in the area. Such complex geological conditions will not only affect local forest growth and surface runoff but also hinder local economic development and human activities. Hence, significant differences in local complex geological conditions and natural disasters influence the change in driver impacts in the Southwestern Mountain.
Mountains around the Sichuan Basin. As shown in Table 4, the proportion of cropland and the relief degree of the land surface are the main factors responsible for the variation in the regression coefficients among the four driving forces in the Mountains around Sichuan Basin. Cropland is one of the dominant land classes in these regions. It is worth noting that the local farmland is mainly located on steep slopes, especially within the slope range of 15 to 35, as a result of the local topographical environment [80]. This phenomenon is more prevalent in regions such as Central Hill and Northern Hill. Slope farmland can lead to soil nutrient loss, soil structure deterioration, soil erosion, and increased surface runoff, which in turn pose a threat to local vegetation growth and economic and social development. Hence, the potential threat of sloping farmland to local drivers should be taken into account when managing the variations in driver influences in the Mountains around Sichuan Basin.
Low Mountains and Hills. Tables 5 and 6 indicate that the proportion of built-up land, the proportion of cropland, and the relief degree of the land surface contribute to the variations in the effects of drivers on ESs in the Low Mountains and Hills. Since 1991, Sichuan Province has experienced rapid economic development and urban expansion. This phenomenon was prominent in the Low Mountains and Hills, where the velocity of economic development is faster than in other places due to the relatively flat terrain. As a result of economic growth and urbanization, the area has expanded significantly in terms of both built-up land and cropland [33]. During this process, extensive highly artificially modified land may alter the hydrological flow path and soil properties and have long-term consequences on local surface runoff and plant growth [81]. Therefore, significant differences in urbanization levels should be considered when exploring and managing the relationship between ESs and driving forces and promoting ecological sustainability in the Low Mountains and Hills.

Implication
The goal of environmental management is the sustainable utilization of ecosystems to meet human demand [82]. Numerous studies have explored regional management based on the biophysical or economic values of ESs; however, these studies have paid less attention to the regional differences in the effects of driving forces on ESs. According to the results of CCA, the relationship between drivers and ESs varied at different spatial scales, with corresponding shifts in the associated environmental factors within each ecoregion. Therefore, in highly heterogeneous areas, the overall policy is difficult to implement [19]. It is important to fully understand the associations between drivers and ESs across multiple spatial scales and take into account scale-related information to optimize localized ecological management [6].
According to our results, a highly heterogeneous region should be able to incorporate significant differences in topographical conditions (e.g., altitude, slope) when exploring the relationship between drivers and ESs and determining the appropriate environmental management strategy at the provincial level, while differences in physical and socioeconomic conditions (e.g., urbanization level, human activity, local vegetation coverage) should be considered at the ecoregional and subecoregional levels.
More specifically, our results suggest that in Sichuan Province, topographical factors (e.g., the relief degree of the land surface, digital elevation model) should be considered as ecoregion boundaries to facilitate localized natural conservation, forestry management and ecological management. In the Western Sichuan Plateau, the focus of ecological management should be on protecting and recovering local natural vegetation and coordinating the relationship between socioeconomic development and ecological management. In response, a series of ecological restoration projects, such as the "Grain for Green Program (GFGP)" and "Natural Forest Conservation Program (NFCP)", could be implemented in an effort to control and mitigate the degradation of the local environment. In the Southwestern Mountain, managers should prioritize optimizing the coverage and structure of local forests to improve their ecological functions and maintain the sustainability of the local environment. With respect to local complex geological conditions, ecological programs for rocky desertification and soil erosion should be combined to enhance the quality, productivity, and diversity of local ecosystems. In the Mountains around Sichuan Basin, managers should be concerned about the potential threat of local sloping farmland. The implementation of soil erosion control measures (e.g., protection of the topsoil layer, application of organic fertilizer and renovation of sloping farmland) and forest conservation projects (e.g., returning farmland to forest, natural forest protection measures, prohibition of timber harvesting) should be coordinated to enhance the resistance of local ecosystems to water erosion [82]. In the Low Mountains and Hills, local managers should be aware of the negative effects of urbanization and the importance of reasonable urban planning. Improving forest coverage and protecting ecological functional zones are encouraged to optimize the local ecological environment and promote regional sustainable development [83]. In addition, we suggest the need to develop agroforestry to ensure supply services while guaranteeing the integrity and diversity of the region's ecosystems.

Conclusions
Based on analyzing the changes in the effects of driving forces on ESs at three spatial scales, identifying the main factors associated with these variations at each spatial scale, and discussing the implications of localized ecological management by analyzing the relationships between environmental factors and drivers, this study examines the significance of the scale dependency of drivers on ESs for localized ecological management. According to these findings, GDP, POP, runoff, and NDVI all have different effects on soil conservation and water yield at three spatial scales (provincial, ecoregional, and subecoregional scales) in Sichuan Province. The results of the CCA indicated that the major factors associated with the regression coefficients of drivers are significantly different at the three spatial scales (provincial, ecoregional, and subecoregional scales). At the provincial (P) scale, physical, topographical, and biophysical factors were the principal factors associated with the variation in the local regression coefficients of the four driving forces. At the ecoregional and subecoregional scales, physical, socioeconomic, topographical, and biophysical factors all contribute to these variations. Our results suggest that significant differences in topographical conditions (e.g., altitude, slope) can be incorporated for exploring the relationship between drivers and ESs and optimizing ecological management at the provincial scale, whereas significant differences in physical and socioeconomic conditions (e.g., urbanization levels, human activity, local vegetation coverage) are more meaningful for localized management at the ecoregional and subecological scales. These findings highlight the importance of awareness among scientists and decision-makers of scale issues of the spatially varying relationships between driving forces and ESs and provide scale-based references for local policy-makers to optimize ecological management at different scales.